Search for new physics in B to VV decays and other hot topics from Belle 



Katsumi Senyo, for Belle Collaboration 
Department of Physics, Nagoya University, 
Furo-cho, Chikusa-ku, Nagoya, Japan 



We report studies in polarization in B decay into two vector mesons with data equivalent 
to 2531b" 1 on T(45) resonance at KEKB. In B° -> (j)K*°, B+ -» <j>K* + and B+ -> p + K*° 
decays, we determine f L = 0.45±0.05(stat.) ±0.02(syst.), ft = 0.52 ± 0.08(stat.) ± 0.03(syst.) 
and ft = 0.43 ± 0.11(stat.)ig'o!(syst.) respectively, where ft is a ratio between longitudinal 
and transverse polarization, and a naive theoretical estimation assumes fn ~ 1. The discrep- 
ancy from 1 in fh may suggest existence of new amplitude within or/and beyond the Standard 
Model. 



1 Introduction 

Naive factorization in the Standard Model (SM) predicts that the longitudinal polarization 
fraction (f£) in B meson decays to light vector-vector(VV) final states is close to unity^ . In the 
tree dominated B + — ► p + p° and B° — ► p + p~ decays, this prediction has been confirmedEEEEl 
In the contrast, for the pure b — > s penguin B — > (f>K* decay, Bellc^and BABAiPhave found that 
longitudinal and transverse polarization fraction are comparable, which is in disagreement with 
the factorization expectation. Possible explanations for this discrepancy include enhanced non- 
factorizable c ontr ibutions such as penguin annihilatioJP, large SU(3) breaking in form factor^, 
or new physicPM It is therefore important to perform polarization measurements in B — > cfrK* 
with larger data set and in other VV modes, in particular, in the pure penguin b — ► sdd decay 
B + -> p + K*°. 

In this study, we use a amount of 253fb _1 of data on T(4 S) resonance, equivale nt t o approx- 
imately 274.5 x 10 6 BB pairs, collected by Belle detectoffllat KEKB e+e" collideM Detailed 
description of Belle detector is found in elsewherP^. 



2 B -> 4>K* Analysis 

The event reconstruction is performed as follows; Candidate B mesons are reconstructed from 
(f) and K* candidates and are identified by the energy difference AE = Ej^ s — E^ m , the 



beam constrained mass M\, c = J {E^ m ) 2 — (p^f 18 ) 2 , and K + K invariant mass (M K + K -). 
^beam * s ^ ne beam energy in the center-of-mass (cms), and E^ s and p c ^ s are the cms energy 
and momentum of the reconstructed B candidate. The i?-meson signal region is defined as 
M bc > 5.27GeV/c 2 , \AE\ < 45MeV, and \M K+K - - M^\ < 10MeV/c 2 . The invariant mass of 
the K* candidate is required to be less than 70MeV/c 2 from the nominal K* mass. The signal 
region is enlarged to — lOOMeV, AE < 80MeV for B+ -> 4>K* + (K* + -> A+vr ) because of the 
impact of shower leakage on AE resolution. The additional requirement cos 6k* < 0.8 is applied 
to reduce low momentum it background, where Ok* is the angle between the direction opposite 
to the B and the daughter kaon in the rest frame of A*. 

The dominant background is e + e~ — > qq{q = u, d, c, s) continuum production. Several 
variables including .S'J-^T the thrust angle, and the modified Fox- Wolfram moments^^ are used 
to exploit the differences between the event shapes for continuum qq production (jet-like) and for 
B decay (spherical) in the cms frame of the T(4S). These variables are combined into a single 
likelihood ratio 1Z S = C S /(C S + C q q, where C s (C g g) denotes the signal ( continuum) likelihood. 
The selection requirements on 1Z S are determined by maximizing the value of N S /^N S + in 
each i?-flavor-tagging quality regioJt^l, where N s (Nf } ) represents the expected number of signal 
(background) events in the signal region. 

Backgrounds from other B decay modes such as B — > K + K~K*, B — > /o(980)A*(/o — > 
A + A~), B — > (fiKir, B — * K + K — Kir, and cross-feed between the <pK* and cpK decay channels 
are studied by Monte Carlo(MC) simulations. The uncertainty in the /o(980) width (40-100 
MeV/c 2 pID i s ta k en 

as a source of systematic error. The contribution from B — > K + K K* 
(B — > foK*) is estimated to be 1-7% (1-3%) of the signal yield. The background from B — > cfiKir 
decays is evaluated with fits to the Kir invariant mass and is found to be about 1%. To remove 
the contamination from <f>K decays, these decays are explicitly reconstructed and rejected. 

The signal yields (Ns) are extracted by extended unbinned maximum-likelihood fits per- 
formed simultaneously to the AE, M^ c and M k +k~ distributions. Reconstructed B candidates 
with \AE\ < 0.25GeV, M bc > 5.2GeV/c 2 , and M K + K - < 1.07GeV/c 2 are included in the fits. 
The signal probability density functions (PDFs) are products of Gaussians in AE and M^c, 
and a Breit-Wigner shape in M k +k~- Bifurcated Gaussians (Gaussians with different widths on 
either side of the mean) are added to model the tails in the AE distributions. The means and 
widths of AE and M& c are verified using B — > J/tpK* decays. The mean and width of the <j> 
mass peak are determined using an inclusive eft — > K + K~ data sample. 

The PDF shapes for the continuum events are parameterized by and ARGUS functioJt^ in 
M^c, a linear function in AE, and a sum of a threshold function and a Breit-Wigner function in 
Mk+k~- The parameters of the functions are determined by a fit to the events in the sideband. 
The signal and background yields are allowed to float in the fit while other PDF parameters are 
fixed. The direct CP asymmetries, Acp = Z^)^jv(i?^>/) ' are a ^ so s t uc bed. The measured 
signal yields and direct CP asymmetries are summarized in Tabled The distributions of AE, 
Mi, c and M k +k~ are shown in Figure ^ 

The decay a ngle s of a iJ-meson decaying to two vector mesons <f> and K* are defined in the 
transversity basis^. The x — y plane is defined to be the decay plane of K* and the x axis is in 
the direction of the cj> meson. The y axis is perpendicular to the x axis in the decay plane and 
is on the same side as the kaon from the K* decay. The z axis is perpendicular to the x — y 
plane according to the right-hand rule, #tr(<^tr) is the polar(azimuthal) angle with respect to the 
z-axis of the K + from <f> decay in the (j) rest frame, and 9k* is defined earlier. 
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Figure 1: Distributions of the AE, M^ c and 
M k +k- for B° -> 0A"*°((a), (b) and (c)), and for 
B + -* 4>K* + ((d), (e) and (f)), with other variables 
in signal region. Solid curves show the fit results. 
The continuum background components are shown 
by the dashed curves. The dark shaded areas rep- 
resent the contributions from B — > K + K~ K* and 
B — > foK* decays. 



Figure 2: Projected distributions of the three 
transversity angles for B° — > <pK*°((a), (b) and 
(c)), and for B + -> <j}K*+ ((d), (e) and (f)). Solid 
lines show the fit results. The points with error 
bars show the efficiency corrected data after back- 
ground subtraction. The two K*° decay modes 
are combined in (d), (e) and (f). The discontinu- 
ity in (d) is due to the requirement of cosOk* in 



B + -> <f>K* + {K* 



K 



Table 1: Number of events observed in the signal region(JV e „) , signal yields (N s ), and the direct CP 
asymmetries(j4cp) obtained in the fits, with statistical and systematic uncertainties. 



Mode 


N ev 


N s 


A C p 




309 


173 ± 16 


0.02 ±0.09 ±0.02 


<f>K* + 
<pK*+ 


(K+ir°) 


173 
76 
97 


37.9±?; 7 
47.3±|;i 


-0.02 ±0.14 ±0.03 
-0.14 ±0.21 ±0.04 
0.09 ±0.19 ±0.04 



Table 2: The decay amplitudes obtained for B° — * <f>K*° and B + — > 4>K* + . the first uncertainties are statistical 

and the second are systematic. 



N ev 


0K* U <f>K*+ 


\Ai_\ 2 
arg(A||)(rad) 
arg(j4^)(rad) 


0.45 ± 0.05 ± 0.02 0.52 ±0.08 ±0.03 
0.30 ± 0.06 ± 0.02 0.19 ±0.08 ±0.02 
2.39 ±0.24 ±0.04 2.10 ± 0.28 ± 0.04 
2.51 ±0.23 ±0.04 2.31 ±0.30 ±0.07 



Table 3: A and E values obtained from the decay amplitudes measured for B and B° separately. 





B° 


B0 


I-Aol 2 


0.39 ± 0.08 ± 0.03 


0.51 ±0.07 ±0.02 


\A ± \ 2 


0.37 ±0.09 ±0.02 


0.25 ±0.07 ±0.01 


arg(A||)(rad) 


2.72±g;^ ±0.14 


2.08 ±0.31 ±0.04 


arg(Aj_)(rad) 


2.81 ±0.36 ±0.11 


2.22 ±0.35 ±0.05 




0.13+^4 ± 0.04 


0.28 ±0.08 ±0.01 




0.03 ±0.08 ±0.01 


0.03 ±0.06 ±0.01 



Table 4: The measured decay amplitudes and triple-product correlations in the B° and B° samples. 



Aoo = 


0.45 ±0.05 ±0.02 


^oo — 


-0.06 ±0.05 ±0.01 


A llll = 


0.24 ±0.06 ±0.02 


S|||| = 


-0.01 ±0.06 ±0.01 


A_ll = 


0.31 ±0.06 ±0.01 


E±± = 


-0.06 ±0.05 ±0.01 


S±o = 


-(Ulio;" ± 0.04 


A± 


= 0.16±g;J2 ± 0.03 


= 


-0.06 ±0.10 ±0.01 


Ax)| = 


= 0.01 ±0.10 ±0.02 


A|| = 


-0.45 ±0.11 ±0.01 


%> = 


-0.11 ±0.11 ±0.02 
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where Aq, An, and A± are the complex amplitudes of the three helicity states in the transversity 
basis with the normalization condition |Ao| 2 + |^4||| 2 + |^._l| 2 = 1> an d t] == 1(— 1) corresponds 
to B(B) mesons and is determined from the charge of the kaon or pion in the K* decay. The 
longitudinal polarization component is denoted by Aq and Aj_(Aii) is the transverse polarization 
along the z-axis (y-axis). The value of \A±\ (|A)| 2 + l^||| 2 ) is the CP-odd (CP-even) fraction 
in the decay B — > (j)K*°^i. The presence of final state interactions (FSI) results in phases that 
differ from either or ±tt. 

The complex amplitudes are determined by performing an unbinned maximum likelihood fit 
to the B — ► <f)K * candidates in the signal region. The combined likelihood is given by 



£ = II e (^ tr ' COS ^ tr ' cos ® K *) X] /j-^jC^tr, cos tr , cos 8 K * 



(1) 



where j denotes the contributions from 4>K*, qq, K + K K* and foK*; Rj is the angular distri- 
bution function (ADF). The ADF R q q is determined from side-band data, and Rk+k-k* from 
events with 1.04GeV/c 2 < M K+K - < 1.075GeV/c 2 ; R foK * is obtained from B f K* MC. 
The detection efficiency (e) is determined using MC simulations assuming a phase space decay. 
The fractions are parameterized as a function of AE and M k +k~- The value of arg(Ao) 
is set to zero and \An\ is calculated from the normalization condition. The four parameters 
(l-Aol, \A±\, arg(j4||), and arg(j4^)) are determined from the fit. There is a two-fold ambiguity in 

the solutions for the phases; the chosen set of solutions is the one suggested in the reference!^! 
Figure El shows the angular distributions with the projections of the fit superimposed. The 
obtained amplitudes are summarized in Tabic [21 

The systematic uncertainties on the amplitudes are dominated by the efficiency modeling 
(4-5%), continuum background (3-4%), slow pion efficiency (2-3%), and K + K~K* ADF (1- 
2%). The remaining possible systematic errors such as the angular resolution, signal yields, 
background from higher K* states, and width of the /o are estimated to be less than 1%. 

The triple-product for a B meson decay to two vector mesons takes the form q ■ (e| x ej), 
where q is the momentum of one of the vector mesons , and e{ and €2 are the polarizations of 
the two vector mesons. The following two T-o quantities 



A° 



Iui(A ± Aq), A\, = lm(A ± AZ 



(2) 



provide information on the asymmetry of the triple products. The SM predicts very small values 
for Aj, and Aj,. The comparison of these triple product asymmetries (Aji and Aj,) with the 



corresponding quantities for the CP-conjugated decays (Aj, and Ay) provides an observable 
sensitive to T-violation. 

Additional variables that can be measured through angular analysis are sugg 

esteJpl and are 

given by 

A ±i = -Im(A ± A* - A±A*), A,| = Re(A||A£ + AyA^), 
= -Im(A ± A* + A X A*), E| j0 = Re(A,|A^ - AyA^), 



Aaa = t,(|^a| 2 + |A A | 2 ), £aa = ±(|A A | 2 - |A A | 2 ) 



(3) 



where the subscript A is one of 0, ||, or _L and i is one of or ||. The variables Aj_o and Aj_|| are 
sensitive to T-violating new physics. The following equations should hold in the absence of NP: 

S AA = 0,S|| = 0,A ±i = 0. (4) 

By separating B° and B° samples and rearranging fitting parameters in the unbinned max- 
imum likelihood fit, we obtain the decay amplitudes for the B° and B°, the triple-product 
correlations, and the other NP-sensitive observables, which are given in Table 01 and 0J 

In summary improved measurements of the decay amplitudes B — > <f>K* are presented, 
based on fits to angular distributions in the transversity basis. The results are consistent with 
our previous measurements^ but with improved precision. The measured value of |Aj_| 2 shows 
that CP-odd (|Aj_| 2 ) and CP-even(| Aq\ 2 + |A||| 2 ) components are present in (j)K* decays in a 
ratio of 1:2. Both phases of Aj_ and An differ from zero or — tt rad by 4:3 standard deviations 
(cr) which provides evidence for the presence of final state interactions. The measured direct 
CP asymmetries in these modes are consistent with zero. These correspond to 90% confidence 
level limits of -0.14 < A C p{4)K* (K + 'K~)) < 0.17, and -0.25 < A C p{4>K* + ) < 0.22. Difference 
between triple products asymmetries (Aj_ 0j |j) which are sensitive to T-violation are consistent 
with zero. The equations S AA = 0, £|| = 0, and Aj_j = should be hold in the absence of 
NP. Our data does not show any significant violation of these relations. Measurements of the 
T-violation sensitive differences between triple product asymmetries, A T — A r and Aj, — A T , 
indicate no significant deviations from zero. Our data indicates no significant deviations from 
the expectations: S AA = 0, Smq = 0, and Aj_j = 0, indicating no evidence for new physics. 



3 B -> pK* Analysis 

We select B + — > p + K*° candidate events by combining three charged tracks (two oppositely 
charged pions and on kaon) and one neutral pion. Each charged track is required to have 
a transverse momentum px > O.lGeV/c and be originated from the interaction point (IP). 
Candidates n° mesons are reconstructed from pairs of photons that have an invariant mass in 
the range 0.1178GeV/c 2 < M^o < 0.1502GeV/c 2 . The tt° candidates are kinetically constrained 
to the nominal 7r° mass. In order to reduce the combinatorial background, we only accept ir° 
candidates with momenta p n o > 0.40GeV/c in the cms. Candidate p + mesons are reconstructed 
via their p + — > 7r + 7r° decay, and the 7r + 7r° pairs are required to have an invariant mass in 
the region 0.62GeV/c 2 < M^+^o < 0.92GeV/c 2 . Candidate K*° mesons are selected from the 
K*° -> K + tt~ decay with an invariant mass 0.83GeV/c 2 < M K+7T - < 0.97GeV/c 2 . To isolate 
the signal we accept events in the region M bc > 5.2GeV/c 2 and -0.3GeV < AE < 0.3GeV, and 
define a signal region in M bc and AE as M bc > 5.27GeV/c 2 and -O.lOGeV < AE < 0.06GeV 
respectively. The continuum process is the main source of background to be suppressed. In 
addition to the continuum process reduction at B — > (j)K* analysis, the displacement along the 
beam direction between the signal B vertex and that of the other B, 5z, also provides separation. 




M bc (GeV/c 2 ) AE (GeV) 



Figure 3: Projections of Mbc for events in the AE 
signal region (left), and projection of AE in the 
A'hc signal region (right). The solid curves show the 
results of the fit. The hatched histograms represent 
the continuum background. The sum of the b — > c 
and continuum background component is shown as 
dot-dashed lines. 
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Figure 4: Signal yields obtained from the Mbc — AE 
distribution in bins of M^+ n o (left) for events in the 
K*° region and in bins of Mk-k (right) for events in 
the p region. Solid curves show the results of the 
fit. Hatched histograms are for the non-resonant 
component. 



For B events, the average value of 5z is approximately 200/im while continuum events have a 
common vertex. This suppression removes 99.3% of the continuum background while retaining 
41 % of the B + — > p + K*° events. The MC-determined efficiency with all selection criteria 
imposed is 2.7% for longitudinal polarization(Ao) and 4.0% for transverse polarization (^4±). 

To investigate backgrounds from b — > c decays, we use a sample of BB MC events correspond- 
ing to an integrated luminosity of 412fb _1 . We find a contribution from B + — > D® (K + tt~ it q )tt + 
decays in the p or K* sideband region and require) M Kn7T o - M D o\ > 0.050GeV/c 2 to veto 
these events. Among the charmless B decays, potential background arise from B + — > a®K + , 
B + — ► p + K*°(1430), non-resonant B + — ► p + K + n~ and B + — ► iv'*°7r + 7r°. We separate signal 
from these backgrounds by fitting the p and K* invariant mass distributions. 

We extract the signal yield by applying an extended unbinned maximum-likelihood fit to the 
two-dimensional M^ — AE distribution. The fit includes components for signal plus backgrounds 
from continuum events and b — > c decays. The PDFs for signal and b — > c decay are modeled 
by smoothed two-dimensional histograms obtained from large MC samples. The signal PDF is 
adjusted to account for small differences observed between data and MC for a high-statistics 
mode containing ir° mesons, B + — ► D (i\~ + 7r~7r )-7r + . The continuum PDF is described by a 
product of a threshold (ARGUS) function for Mb c and a first-order polynomial for AE, with 
shape parameters allowed to vary. All normalizations are allowed to float. Figure 01 shows the 
final event sample and the fit results. The five-parameter (three normalizations plus two shape 
parameters for continuum) fit yields 134.8 ± 16.9 B + — * K + 7r~7r + 7T° events. 

We further distinguish the p + K*° signal from non-resonant decays such as B + — > p + K + -K~ 
or B + — ► K*°7r + 7r° by fitting the M n + n o and Mk-k invariant mass distributions. The signal yields 
obtained from the M^ — AE fit for different M^+^o and Mk-k bins are plotted in Figure^ where 
the M^+^o distribution is for events in the K*° region (0.83GeV/c 2 < M Kn < 0.97GeV/c 2 ) and 
the the Mk-k distribution is for events in the p region (0.62GeV/c 2 < M^+^o < 0.92GeV/c 2 ). 
We perform separate x 2 fits to the M wn or Mk-k distributions. Each fit includes components for 
signal and non-resonant background, the signal p and K* PDFs are modeled bv relativistic P- 
wave Breit-Wigner functions with means and widths fixed at their known values'*"", the PDFs are 
convolved with a Gaussian of a = 5.3MeV, which is obtained by fitting the D°(K~ir + ) invariant 
mass, to account for the detector resolution. The non-resonant component is represented by 
a threshold function with parameters determined from MC events where the final states are 
distributed uniformly over phase space. The M nn mass fit gives 125.4 db 15.8/0 and —0.3 db 3.0 
non-resonant events in the K* mass region. The statistical significance of the signal, defined as 
\Jxo~ Xmin' wri ere Xmin i s the x 2 value at the best-fit signal yield and Xo ^ s the value with the 
K*° signal yield set to zero, is 5.3o"(5.2cr with the inclusion of systematics). The contribution 
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Figure 5: Fit to background-subtracted helicity distributions. The solid histograms show the results. The dot- 
dashed (dashed) histograms are the Ao(A±) component of the fit; the dotted histograms are for non-resonant 
pKn. The low event near cos(9 p = 1 is due to the p^o > 0.4GeV/c requirement. 



from non-resonant p + K + ir~ is significant and is taken into account in both the branching fraction 
and polarization determinations, while we neglect the non-resonant K*°tt + tt° contribution. 

We use the p + — ► 7r + 7r and K*° — > K + ir~ helicity-angle (0 P , Ok*) distributions to determine 
the relative strengths of | viol 2 an d |^4±| 2 - Here p (0k*) is the angle between an axis anti-parallel 
to the B flight direction and the ir + (K + ) flight direction in the p(K*) rest frame. For the 
longitudinal polarization case, the distribution is proportional to cos 2 9„ cos 2 Ok* , and for the 
transverse polarization case, it is proportional to sin 2 9 p sm 2 K * P Figure Q shows the signal 
yields obtained from M^ c — AE fits in bins of the cosine of the helicity angle for p and K* . We 
perform a binned simultaneous x 2 fit to the p and K* helicity-angle distributions. The fit includes 
components for signal and non-resonant pKix. PDFs for signal Aq and A± helicity states are 
determined from the MC simulation. The helicity-angle distribution for data in the high Mk-k 
sideband region l.lGeV/c 2 < Mk-k < 1.5GeV/c 2 , where pKir events dominate, is consistent 
with a cos 2 0-like cos P and a flat cos Ok* distribution. Thus, we assume an S-wave Kir system 
and model the non-resonant B — > pKn PDF based on the MC simulation. The fraction of the 
non-resonant component is fixed at the values obtained from the K* mass fit. The two parameter 
(normalizations for Aq and A±) fit result deviates from 100% longitudinal polarization with a 
significance of 4.9a (4.3a including systematic uncertainties). The significance is defined as 

\j Xq ~ Xmin' where Xmin * s the x 2 value at the best-fit and Xo * s the value with the longitudinal 
polarization fraction set to 100%. 

The largest uncertainties in the polarization measurement are due to uncertainties in the non- 
resonant pKir PDF, potential scalar-pseudoscalar (S — P) interference, and the non-resonant 
fraction. We assign 1q 0,3 % systematic error for the non-resonant PDF. This uncertainty is 
estimated by adding a 1/3 flat component to the p helicity PDF for non-resonant pKir in 
the helicity fit. Interference of the longitudinal amplitude Aq with the <S-wave (Kir) system 
introduces a term with a 2e tA ^\A p Kn\ cos Ok* dependence, where Acf> is the phase difference and 
|^4pif7r| is amplitude of the B — > pKir decay. The S — P wave interference disappears in the cos P 
distribution, which is integrated over cos Ok* ', however it remains in the cos Ok* distribution. We 
include an additional linear function for the interference term in the cos Ok* helicity and redo 
the x 2 fit. The resulting small change in fi, 0.5%, is assigned as the systematic uncertainty for 
the S — P interference, systematic error is assigned for the uncertainty in the fraction of 

non-resonant pKir, obtained by varying the non-resonant fraction by ztlu. Adding the various 
systematic error contributions in quadrature, we obtain the longitudinal polarization fraction in 
B+ -> p+K*° decays, 

f L (B+ - p + K*°) = 0.43 ± 0.11(stat.)lH|(syst.). (5) 

To calculate the B + — > p + K*° branching fraction, we use the invariant mass fit result and 
MC-determined efficiencies weighted by the measured polarization components. We consider 



systematic errors in the branching fraction that are caused by uncertainties in the efficiencies 
of track finding, particle identification, 7r° reconstruction, continuum suppression, fitting, polar- 
ization fraction. We assign an error of 1.1% per track for the uncertainty in the track efficiency. 
This uncertainty is obtained from a study of partially reconstructed D* decays. We also assign 
an uncertainty of 0.7% per track on the particle identification efficiency, based on a study of 
kinematically selected D* + — ► D°(K~tt + )tt + decay. A 4.0% systematic error for the uncertainty 
in the 7r° detection efficiency is determined from data-MC comparisons of rj — > 7r 7r 7r° with 
r\ — ► 7T + 7r~7r° and n — > 77. A 4.5% systematic error for continuum suppression is estimated 
from studying the process B + — > D (K + ir~TT )ir + . A -4.2%/+1.7% error due to the uncer- 
tainty in the fraction of longitudinal polarization is obtained by varying fi by its errors. The 
uncertainty in non-resonant K*mr background gives a contribution of -2.2%/+0% in addition 
to -3.0%/+2.3% error from uncertainties in the background from other rare B decays. A 1.1% 
error for the uncertainty in the number of BB events in the data sample. A 7.1% error for 
possible bias in the x 2 fit^l obtained from a MC study is also included. The quadratic sum of 
all of these errors is taken as the total systematic error. We obtain the branching fraction 

B{B + -> p + K*°) = (8.9 ± 1.7(stat) ± 1.2(syst) x 10~ 6 . (6) 

In summary, we have observed the B + — ► p + K*° decay with a statistical significance of 5.3a. 
We measure the branching fraction to be (8.9 db 1.7(stat) ± 1.2(syst) x 10~ 6 . We also perform a 
helicity analysis and find a substantial transversely polarized fraction with a statistical signifi- 
cance of 4.9cr. The longitudinal polarization fraction fi, measured is similar to the surprisingly 
low value found in b — > sss decays B — ► cf)K * . 
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